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ABSTRACT 

Filled arrays of bolometers are currently being employed for use in astronomy from 
the far-infrared through millimeter parts of the electromagnetic spectrum. Because of 
the large range of wavelengths for which such detectors are applicable, the number of 
modes supported by a pixel will vary according to the specific application of a given 
available technology. We study the dependence of image fidelity and induced polariza- 
tion on the size of the pixel by employing a formalism in which diffraction due to the 
pixel boundary is treated by propagating the second-order statistical correlations of the 
radiation field through a model optical system. We construct polarized beam pattern 
images of square pixels for various ratios of p/X where p is the pixel size and A is the 
wavelength of the radiation under consideration. For the limit in which few modes are 
supported by the pixel (p/X < 1), we find that the diffraction due to the pixel edges is 
non-negligible and hence must be considered along with the telescope diffraction pattern 
in modeling the ultimate spatial resolution of an imaging system. For the case in which 
the pixel is over-moded (p/X > 1), the geometric limit is approached as expected. This 
technique gives a quantitative approach to optimize the imaging properties of arrays of 
planar detectors in the few- mode limit. 

Subject headings: techniques: polarimetric — techniques: photometric — methods: 
statistical 



1. Introduction 



Filled arrays of planar detectors have a long history of astronomical use in the optical (CCDs) 
and infrared (InSb, HgCdTel) regions of the electromagnetic spectrum. In recent years, planar 
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bolometers have enabled the extension of this technique into the millimeter range where coherent 
detection techniques have traditionally been used. Examples of planar millimeter-wave elements 
can include both patches of absorbing material, which can be multimode if made large enough, 
or single- mode antenna-coupled structures. In this paper, we are concerned with detectors that 
operate through the absorption of radiation in a lossy thin-film or bulk material. 

The constituent detectors of planar arrays are normally considered in the geometric limit. In 
other words, the coherence length associated with the loss mechanism that is responsible for power 
absorption is small compared with the detector's physical size. In this limit, the only diffractive 
effects considered are those induced by the optical system that controls the illumination of the 
detectors, and the individual pixels of the array merely extract power from the incoming field 
according to their geometrical size. In such a model, the sampling of a telescope-induced Airy disk 
is assumed to be perfect, and the Airy disk sets the angular resolution limit of the optical system 
as long as there are a sufficient number of pixels within the disk. 

In order for the detector to behave in a geometrical manner, the ratio of the pixel size (p) 
to the wavelength (A) must be significantly greater than unity. In many cases, however, the indi- 
vidual pixels are nearly a wavelength in size, and so their behavior cannot be described in terms 
of geometrical shape alone, as the pixel boundary conditions in this limit become important. For 
example, the currents flowing around the edges of the individual pixels will lead to beam patterns 
that exhibit diffractive phenomena, including potentially complicated states of polarization. An- 
other way of describing the change from geometrical to diffractive behavior is that the number of 
electromagnetic modes that can be supported by a detector decreases as the size of the detector 
decreases. 

In millimeter-wave astronomy, practical limitations such as the size of the optical system and 
those associated with physical properties relating to the speed of bolometric detectors (thermal 
conductivity and heat capacity), motivate the construction of fast optical systems that have lower 
p/X than their high-frequency counterparts. In this case, the diffractive effects of the pixels and their 
associated architecture are no longer negligible, and the consequences for the imaging capability of 
the sys tem must be und erstood. Such image degradation is analogous to that of infrared photodiode 
arrays (jHollowavl Il986l ) ; however, in this case, since 1 < p/X < 5, the image degradation is not due 
to diffractive effects, but rather to diffusive spreading of charge carriers in a common semiconductor 
absorber. In both cases, these focal plane effects must be considered when designing the front-end 
optics. A partial list of current and planned instruments that e mploy planar bolom eters is given 
in Table ffl These includ e HAWC /SOFIA dHaroer et all l2004h, SHARC II/CSO (bcwejletal 



200 3), SCUBA 2/JC MT (Holland et al 



( Staguhn et alJl200fl ). and ACT (IFowler 



2004 ) 



2009), PAR/GBT (IDicker et alJl200fl ). GISMO/IRAM 



In this work, we explore the diffractive effects of pixels for which 0.25 < p/X < 4.00. We adapt 
a method in which the statistical correlations of an incoming radiation field, a s distinct from the 



field itself, are propagated through a model optical system (jWithington 



2006 



Withington et al 
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2003). Care has been taken in the simulations to ensure that the diffractive effects seen on the sky 
are dominated by the pixels and not by the optical system that illuminates the array. This has 
been done with the specific purpose of isolating the effects of using small pixels from the effects of 
the telescope optics. We find that when p/A > 1, the geometric limit is recovered, but that when 
p/A < 1, diffractive effects are non-negligible, and must be taken into account when designing 
instruments that employ such pixels. 



In Section 2, we briefly review the method of IWithington et al.l (|2003l ) as it applies to our 
calculations. In Section 3, we describe the results of the calculated polarized and unpolarized beam 
patterns for various values of p/A. In Section 4, we suggest a method for using the beam patterns to 
calculate the cross-coupling between pairs of pixels when viewing an extended, incoherent source. 
In Section 5, we briefly consider the beam shape of a polarized pixel. Finally, conclusions and 
future work are discussed in Section 6. 



The Method 



Sections 2.1 and 2.2 are provided as a brief summary of IWithington et al.l (|2003l ). For a more 
thorough discussion of the method they have developed, the reader is invited to consult their work. 
In Section 2.3 we connect the correlation dyadic and the Stokes parameters. We later employ 
this connection in our analysis and presentation of the diffractive effects of the pixels. Finally, 
Section 2.4 concludes with a discussion of the number of modes required in order to ensure that the 
diffraction in the model system is dominated by the pixels and not due to other optical elements 
in the system. 



2.1. Basic Formalism 

The correlations between the components of an electric field at two locations (n and F2) are 
conveniently expressed in terms of a space-domain correlation dyadic: 



E(rx,f 2 ) = {E{r 2 )E (ri)^. (1) 

The angle brackets denote an ensemble average. 

The space-domain correlation dyadic encodes the second-order statistics of the radiation field. 
For purposes of this work, we take the z-axis to be the axis of the optics. The second-order statistics 
can be propagated from one plane (z = constant) to another in the optical system by implementing 
the appropriate boundary conditions. 



Following IWithington et al.l (|2003l ). the electric field can be decomposed into a complete, or- 



thogonal set of vector plane- wave fields, 



*Pi(h,r t ) = ^{k t )e^e^ z e-^\ (2) 
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where, k = kk = k x x + k y y + k z z = kt + k z z, and ii(kt) is the set of unit polarization vectors in 
which i^{kt) = k, the direction normal to the planes of constant phase. Figure Q] illustrates the 
relationships between the relevant vectors. 

Given this basis, and denoting the Fourier amplitude for the mode characterized by transverse 
wave vector k t for the ith polarization as ai(k t ), the electric field can be expressed in as 

n i=l 

It is convenient, and equivalent, to express the correlations of the electric field in terms of the 
fc-domain correlation dyadic, 

_ 3 3 

A(k' t ,k t ) = ^ ^ (°i(*t)aj$)) ?i(h)£M)- (4) 

t=l j=l 

Because we are modeling detectors that have a planar geometry, we are concerned with correlations 
between points in the same plane. That is, z\ = Z2 = 0. Defining fa = x\x+yiy and ft2 = X2X+y2y, 
the space- and fc-domain correlation dyadics are related by two 2-dimensional Fourier transforms. 

= A (%,-k t ) = ^ / f^n^^-^^n^n (5) 



2.2. Model Optical System 

The above formalism is useful for analyzing diffraction in a generic optical system. A diagram 
of the system used to study the behavior of planar bolometer arrays is shown in Figure [21 which is 
a modified version of Figure 1 of Withington et al. (2003). In this model, a pixel is represented by 
a square aperture that is illuminated from the back by a blackbody source of infinite spatial extent. 
The aperture scatters the radiation from the source, thereby introducing additional correlations into 
the reception pattern. The aperture plane is then imaged at the reimaging screen. The reimaging 
screen in this setup represents the projection of the pixel on the sky. 

A key assumption is that it is appropriate to model a planar absorber as a hole in front of a 
blackbody emitter. It is our belief, that for the purposes of understanding general behavior, this 
assumption is quite valid, since to the extent that the absorber is black, the behavior of the system 
is insensitive to the physical realization of the detector. Various small changes can be made to 
take into account the details of a specific absorber, but a generic model that makes no assumptions 
about the precise construction of the pixels, but which brings out key aspects of behavior, is highly 
beneficial. Future simulations can be much more sophisticated. 
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Because the size of the whole imaging array is finite in practice, we work with a discrete basis 



set, 



L 



(7) 



Here, L is the size of the blackbody emitter (assumed to be square), and m and n are the discrete 
rs in tl 

k 2 - h2 < 



mode numbers in the x and y directions, respectively. In addition, k x (m) = k y (n) = ^jP-, and 



k 2 (m,n) = k 2 — k 2 (m) — k 2 (n). 

In this basis, the s pace-domain correlation dyadic for a finite blackbody source is given by 
Withington et al.1 kooi ) as 



E rs (n,f 2 ) = -L Ao (k)J2 (k26rS k2 krks K ^ m ^-^e^ 

m,n 



(8) 



Here, the indices r and s correspond to the components of the space-domain correlation dyadic 
when projected into the (x,y,z) coordinate system. Also, A Q (k) is the amplitude of the electric 
field having wavenumber k and 5 rs is the Kronecker delta. Th ese spatial correlatio ns are the 
combination of those inherent to p ropagating blackbody r adiation (jMehta fc Wolall964l ) and those 
due to the sampling of the source ([Withington et al.l l2003). 



At this point, the /c-domain correlation dyadic can be calculated using Equation [SJ The effect 
of the pixel aperture on the /c-domain correlation dyadic can be expressed in matrix notation as 
A' = SAS^. The /c-domain correlation dyadic to the left of the pixel aperture is A and that to the 
right is A'. The scattering matrix appears twice in the transformation because A consists of two 
electric fields, both of which are scattered by the aperture. 



The scattering matrix elements are given by 

-2 



Sij(r, s;m,n) 



L 2 



ii(m,n) ■ ij(r, s) j (m 



Ml 



(n-sfj 



(9) 



Here, j (%) = sinx/x. This scattering matrix describes the mapping of an incident mode char- 
acterized by (m,n) to scattered mode (r,s), where i and j represent the scattered and incident 
polarizations, respectively. Here, p is the length of a side of a square pixel. 



As discussed in IWithington et al.l ([2003), the pupil stop can be modeled, to first order, simply 
by cutting off the number of modes allowed to propagate through the system. In fact, this simplicity 
is one of the strengths of working with the /c-domain dyadic. However, since we are interested in 
isolating the effect of pixel diffraction, we wish to minimize the diffractive effects of the fore optics. 
Because of this, we ensure that the pupil passes all of the modes produced by the finite-sized 
detector. We then proceed to calculate the space-domain dyadic components using Equation [6] to 
produce an image of the pixel on the reimaging screen. 
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2.3. Image Fidelity and Polarization 

The space-domain correlation dyadic contains all of the information about the second-order 
statistics of the radiation field at any given plane of the optical system. In examining the correlations 
introduced by the optical system at the reimaging screen (or any other plane in the optical system 
for which the space-domain correlation dyadic has been calculated), it is useful to separate the 
correlations into those in which ri 7^ T2 and those in which r = r 1 =¥2- The former describes the 
correlation between the fields in two different locations in the plane that will manifest itself in the 
form of image fidelity. The latter describes the polarization of the beam. 

The polarization can be effectively treated by noticing that the Cartesian components of the 
correlation dyadic can be equivalently expressed in terms of the Stokes parameters. 

I = E xx (r,r) + E yy (r,r) (10) 

Q = E xx (r,r) - E yy (r,r) (11) 

U = m xy (r,r) (12) 

V = QE xy (r,r). (13) 

In this work, we utilize the connection between the space-domain correlation dyadic and the Stokes 
parameters to characterize the diffraction induced by an optical system in terms of its polarizing 
effect on initially unpolarized radiation. We are specifically interested in characterizing both the 
polarization and spatial coherence introduced by the boundary of a pixel. 



2.4. Lateral Size of the System 

The goal of this work is to study the diffractive effects of the pixel edges, and so we would 
like our model optical system to highlight the correlations due to the scattering through the pixel 
aperture, while suppressing those due to other parts of the optics. Because we have disallowed 
mode truncation at the pupil stop, the only way correlations can be introduced into the system by 
anything other than the pixel is through the finite size of the blackbody source, or in other words 
the total size of the region needed numerically to define the Fourier series. We ideally want the 
number of modes to be high enough to approximate a source of infinite lateral size; however, the 
inclusion of more modes slows the calculation. In addition, because the size of the pixel aperture 
will naturally cut off higher modes, inclusion of modes above a certain level is expected to have 
minimal impact on the result. We seek a quantitative way to verify that we have included enough 
modes such that the correlations due to the optics are small compared to those due to the pixel 
aperture. 

To estimate the effect of the fe-domain cutoff introduced by the finite size of the system, we 
ran 1-dimensional calculations for Stokes I (total power) from the center to the edge of the pixel 
for multiple numbers of modes. The results of this study are shown in Figure El In all cases, the 
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number of Fourier modes used in the calculation is significantly greater than the number of modes 
nominally supported by the pixel. Based on these results, we are satisfied that the diffractive beam 
spreading effects that we calculate are due to the electromagnetic boundary conditions at the pixel 
edges, and not from elsewhere in the optical system, or from the finite size of the system used 
to model behavior. The number of modes we have included in our study of pixel diffraction is 
summarized in Table [2j 



3. Characterization of Pixel Diffraction 

We calculate the space-domain correlation dyadic at the reimaging screen for each of five 
different values of p/X. In order to visualize and quantify the effect of pixel diffraction, we then 
calculate the spatial distribution of the Stokes parameters using equations (|10H13p . Both the spatial 
and polarization correlations are captured in this presentation. The calculations of polarized beam 
patterns have the advantage that they provide the capability for a direct comparison to measured 
data. They also can provide guidance when modeling predicted instrument performance. 

Figure H] shows the beam patterns calculated for Stokes parameters I, Q, and U, for various 
choices oip/X. The intensity plots are commonly normalized to the peak flux of the p/\ = 4.00 case. 
Each polarization (Q and U) plot is normalized to the peak flux in the corresponding Stokes I plot, 
and colorbar units are given in per cent. As an aside, we note that Stokes V is also calculable, but 
the symmetry of the system prevents the generation of quadrature correlations between orthogonal 
polarizations in our linear basis. 

In the limit of large p/X, we expect the geometric limit to be recovered. In this case, all of 
the Stokes I flux is contained in the square that defines the physical pixel. Figure U] shows that 
diffractive effects are small at p/X=2, and the geometric limit is indeed nearly recovered by p/X=A. 
In these cases, the diffractive effects of the pixel are more easily seen in the beam patterns of Stokes 
Q and U. For the high p/X cases, Stokes Q traces the regions close to the edges of the pixel in the 
following way: Immediately inside the pixel edge, the polarization direction is parallel to the pixel 
edge. Q is zero at the pixel edge and then switches sign such as to be perpendicular to the pixel just 
outside the edge. Farther out, the polarization is lower and parallel to the pixel edge indicating 
that the highly scattered radiation is polarized parallel to the pixel edge. Stokes U appears at 
the corners of the pixel. Since U is defined as the in-phase correlation between vertical and linear 
polarization, the degree to which Stokes U appears depends on the spatial coherence scale for a 
given wavelength. The point-like behavior may be attributed to current flowing around the corners 
of the detector. As p/X gets large, the area of the beam characterized by non-zero Q and/or U 
becomes increasingly small. We find that ~ 10% of the flux is po larized at p/ X = 1.0, and it drops 



as (p/X) , consistent with the geometrical theory of diffraction ( Keller! 1 1962 ) 



It should be noted that this effect is more severe for U than for Q since for the latter, opposite 
signed regions are located close to one another and are likely to cancel when convolved with a 
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source (or the Airy disk of the telescope). The U beam patterns are spatially separated at high 
p/X, and therefore are more likely to convert unpolarized anisotropics on the sky into polarized 
signals. 

For low p/X, the coherence length is large compared to the size of the pixel. Stokes I transitions 
from a square-like pattern to a circular beam shape. In these cases, the wavelength of the radiation 
is too large to resolve the details of the pixel shape. The effect of diffraction is significant, as 
indicated by the larger spatial extent of Stokes I and the larger values and more extended structure 
of Stokes Q and U. The polarization along the inner edges of the pixel is no longer visible at p/X <1 
as the correlation scale has increased from the multi-mode cases to the single-mode case. 

Figure [5] shows the spatial distribution of fractional polarization (\/Q 2 + U 2 /I) for the extreme 
cases of p/X = 0.25 and p/X = 4.00. In both cases, the on-axis polarization is quite low, and the 
polarization of the scattered radiation quite high. This is expected as scattering and polarization 
are intimately related since both result from the same induced correlations. 



Perhaps the simplest and most important quantity for assessing the effect of p/X on a real 
detector array is the cross-coupling between adjacent pixels. By cross-coupling, we mean the the 
degree to which the emission from one small region of an incoherent sky contributes to the outputs 
of two detectors simultaneously. This is different from the problem of determining the correlations 
between the fluctuations in the outputs of two detectors, which is also possible using the model 
presented, but it is not what we have done here. The loss of resolution introduced by a telescope 
will serve to increase the overlap between the beam patterns of pixels, and thus in some sense, this 
model serves as a best case scenario for given a p/X. In practice, for diffraction-limited systems, 
the imaging system should be designed such that the overlap due to the finite resolution of the 
telescope dominates. 

We will define a cross-coupling factor, which is a function of position on the sky (x, y), as well 
as the two-dimensional separation (Ax, Ay) between the center of the two pixels being compared, 
as 



Here, r = J [I(x', y')] 2 dx' dy', where the integral extends over all space. p(x, y, Ax, Ay) is therefore 
a normalized cross-coupling factor. Figure [6] shows an example of the cross-coupling factor for 
two pixels that are diagonal nearest neighbors. The total cross coupling factor for a uniform sky, 
between two pixels whose beam patterns are separated by Ax and Ay in the x and y directions 
respectively, is simply 



The integral nominally extends over all space, but in reality is limited by the extent of our calculated 



4. Pixel Cross-Coupling 




(14) 




(15) 
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beam profile models. We include a summary of the values of /3(Ax, Ay) for the cases studied 
previously in Fig. [7J 

5. Analysis of Polarized Pixels 

With the interest in astronomical polarimetry growing, a logical extension of planar arrays is 
to pattern the absorbers such that they are sensitive to linear polarization. Doing so allows the 
theoretical possibility of stacking two orthogonally polarized pixels on top of one another in order 
to detect both modes of polarization simultaneously. 

With minor variation (i.e. by eliminating all of the modes of one polarization from the black- 
body source), the analysis technique described here can be used to study the systematic effects one 
would expect when working in the small pixel limit. Figure [8] shows the intensity beam patterns 
for horizontally and vertically polarized detectors (A and B, respectively) for the case of p/X = 0.5. 
The difference between these two images, the incoherent analog of cross-polarization, is shown in 
(C). It is of interest to note that the effective beam size is smeared in the direction of polarization. 
This is consistent with the sign of Stokes Q in Figure |U One might expect that the polarization 
parallel to the pixel edge would be scattered more efficiently; however close to the pixel this is not 
the case. Looking further out in Figure [8] (A and B), one notices that the parallel polarization has 
more support in the highly scattered modes far from the pixel. 

6. Conclusion 

We have explored the consequences of using filled arrays of planar detectors in the limit where 
p/X is small, and found that in this limit, diffraction due to the edge currents in the pixels must be 
considered when designing optical systems. This diffraction has the effect of limiting the angular 
resolution of the instrument for a given plate scale. Thus, by carefully choosing the plate scale, it 
is possible to mitigate this effect. This tends to drive the design such that the pixels oversample 
the Airy function of the telescope. For polarized detectors in this same limit, systematic effects can 
become non- negligible, leading to cross-polarization that is in excess of 10%. 

This model is idealized, but the correspondence of this model to a particular detector imple- 
mentation could be improved. One such improvement involves focusing on the details of how the 
absorbing properties of the pixels are modeled. In our current implementation, the pixels are mod- 
eled as completely incoherent emitting/absorbing apertures. It is possible to tailor this formalism 
to a specific coupling architecture by modeling the detailed material properties of the absorber. In 
addition, for a real array, the pixels near the edge of the array will have different electromagnetic 
boundary conditions than those near the center. This effect could also be included in the model. 

The strength of this technique used in this work is in its ability to handle partially-coherent 
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radiation (jWithington et al.l 120031 ; ICarter fc Somerslll976f). Though coh erent analysis of electro- 



magnetic systems are useful and informative (eg. I Wollack et al.l (|2006l ) and references therein), 
the second-order statistical correlations introduced by these systems are not directly accessible. 
However, consideration of these correlations is essential for understanding and optimizing the per- 
formance of submillimeter and millimeter astronomical systems. 
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Fig. 1. — The coordinate system for the formalism is shown. The (x, y, z) coordinate system is that 
of the optical system with z being the optic axis. The (§i(k t ), ii(k t ), ii(k t )) coordinate system is 
that corresponding to the plane wave mode characterized by kt- The plane of constant phase is 
spanned by £i(fct) and £2(^)5 and the direction of propagation is given by ii(kt) = k. 



-13- 



Table 1. A selection of both current and planned astronomical instruments that employ or will 

employ arrays of planar bolometers 



Instrument 


Array Size 


Detector Type 


A(mm) 


p(mm) 


p/X 


rlAWU/bUr 1 1A 


Izxdz 


Semiconducting Bolometer 


O.Ooo 


1 nn 
1.0U 


18. y 




12x32 


Semiconducting Bolometer 


0.088 


1.00 


11.4 




12x32 


Semiconducting Bolometer 


0.155 


1.00 


6.5 




12x32 


Semiconducting Bolometer 


0.215 


1.00 


4.7 


SHARC II/CSO 


12x32 


Semiconducting Bolometer 


0.350 


1.00 


2.9 




12x32 


Semiconducting Bolometer 


0.450 


1.00 


2.2 




12x32 


Semiconducting Bolometer 


0.850 


1.00 


1.2 


SCUBA 2/JCMT 


64x64 


TES 


0.450 


1.135 


2.5 




32x32 


TES 


0.850 


1.135 


1.3 


PAR/GBT 


8x8 


TES 


3.00 


3.00 


1.0 


GISMO/IRAM 


8x16 


TES 


2.00 


2.00 


1.0 


ACT 


32x32 


TES 


1.13 


1.00 


0.9 




32x32 


TES 


1.33 


1.00 


0.8 




32x32 


TES 


2.07 


1.00 


0.5 



Table 2. A tabulation of the number of modes (along 1 dimension) used for each of the cases 

studied 



p/A 


Number of modes 


0.25 


8 


0.50 


10 


1.00 


12 


2.00 


14 


4.00 


16 
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Fig. 2. — The model optical system that is used for the analysis is shown. The blackbody source 
produces radiation that gets scattered by the pixel aperture. The pupil stop is made large enough 
to pass all of the modes produced by the blackbody, and the radiation pattern is imaged on the 
reimaging screen. 
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Fig. 3. — The effect of varying the number of modes for selected cases of p/A is shown. The y 
coordinate represents the lateral distance from the center of the pixel. The profiles plotted are 
those for total power (Stokes I). 
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Fig. 4. — Stokes I, Q, and U are shown for various values of p/X. Stokes I in each case is normalized 
such that the total power is constant across all five p/X cases such that the common color table. 
Stokes Q and U are normalized to the peak of each of the corresponding Stokes I flux distribution. 
Their corresponding color tables are given in per cent. The pixel size is shown (green squares) for 
the p/X = 0.25 case, and is the same size for all of the cases. 
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Fig. 5. — The fractional polarization (\/ Q 2 + U 2 /I) as a function of position is shown for the cases 
of p/X = 0.25 (A) and p/X = 4.00 (B). The blue square in each image shows the outline of the pixel 
boundary used in each calculation. 
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Fig. 6. — An example plot of the spatial correlation density between two pixels is shown. Here, one 
pixel is displaced relative to the other by 1 pixel unit in each of the x and y directions (p(x, y, 1, 1)), 
and p/X = 1. Contours represent the Stokes I power distributions for the two detectors and are set 
at the following fractions of the peak flux for each pixel's beam: 0.14, 0.28, 0.42, 0.56, 0.70, 0.84, 
and 0.98. 
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Fig. 7. — For the five test cases of p/A, we show the calculated overlap integrals {(3) for pixels with 
the illustrated spatial relationships to the central pixel (for which the overlap integral is defined to 
be 1). 
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Fig. 8. — The beam patterns for horizontally and vertically polarized pixels are shown in A and 
B, respectively. The difference between the two is shown in C, which is analagous to the cross- 
polarization of the pixel. These calculations were done for p/X = 0.5. The size of a pixel is shown 
in the lower left corner of panel A. 



